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We study the small mass limit of the one-loop spinor effective action, comparing the derivative 
expansion approximation with exact numerical results that are obtained from an extension to spinor 
theories of the partial-wave-cutoff method. In this approach one can compute numerically the 
renormalized one-loop effective action, for radially separable gauge field background fields in spinor 
QED. We highlight an important difference between the small mass limit of the derivative expansion 
for spinor and scalar theories. 
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The one- loop effective action in a background field is an important quantity in quantum field theory [H [2]. In 
gauge theories, the one-loop effective action has been calculated exactly and analytically only for certain very special 
gauge field backgrounds, such as those with constant (or, in non-Abelian theories, covariantly constant) field strength, 
based on the seminal work of Heisenberg and Euler [3] , and Schwinger 4J ; and for some very special onc-dimensional 
cases where the field is inhomogeneous [S] . The goal of this paper is to extend further the class of background fields 
for which we can compute the renormalized one-loop effective action. We are particularly interested in studying the 
small fermion mass limit, motivated by its importance for chiral physics, and also by the fact that while the large 
mass limit is well understood 'W , much less is known about the small mass limit for general gauge backgrounds. 

Recently, a new approach has been developed for computing exactly, but numerically, the gauge theory effective 
action in a class of background fields that permit a separation of variables reducing to a set of one-dimensional 
operators. This has been explored in most detail for radially separable backgrounds, where the method has been 
called the "partial- wave cutoff method" [THllJ- The first application of this method was to the full mass dependence 
of the quark determinant in an instanton background, yielding a smooth interpolation between the known large and 
small mass extreme limits [7j. Subsequent papers have concentrated on scalar theories. Indeed, even the instanton 
background computation made use of a special property of self-dual backgrounds that implies that the spinor effective 
action can be expressed directly in terms of the scalar effective action [T^HH]- Here, in this paper, we present the 
spinor approach without assuming self-duality of the background gauge field. 

The basic idea of the "partial- wave cutoff method" is simple: the one- loop effective action requires the logarithm 
of the determinant of an operator, and there is a simple method [known as the Gel'fand-Yaglom method 16]] for 
computing the determinant of an ordinary differential operator, without computing its eigenvalues. For a partial 
differential operator, if the problem is separable down to a set of ordinary differential operators, we can formally sum 
over the Gel'fand-Yaglom results for each term in the separation sum. The technical difficulty is that this sum is 
naively divergent, and this must be addressed by a suitable regularization and renormalization. This has been resolved 
in previous publications for gauge theories with scalars and for self-interacting scalar theories. Here we show how this 
works for spinor theories with background gauge fields that are radially separable, and not necessarily self-dual. 

While this class of radially separable backgrounds is large, and includes important special cases such as instantons, 
monopoles and vortices, there are of course many other physically important background fields that are not separable 
in this way. In these cases we must resort to approximation methods in order to compute the effective action. In 
order to investigate the region of validity of these approximations, we compare our exact numerical results with two 
such approximations, the large mass expansion and the derivative expansion, and show that for spinor theories a new 
feature arises when evaluating the derivative expansion approximation for light fermions. 



We begin with a chiral decomposition of the effective action. The Euclidean one-loop effective action in spinor QED 
is the logarithm of the determinant of the corresponding Dirac operator: 



I. INTRODUCTION 



II. PARTIAL- WAVE DECOMPOSITION OF THE DIRAC OPERATOR 
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Here P = 7^(9;i + iA^{x)) is the Dirac operator in Euclidean 4-diniensional spacetime, and Afj^{x) is the classical 
background gauge field. We will use the following standard representation for the 4x4 Dirac matrices in Euclidean 
spacetime [13]: 



where 



a 



and the ai are the 2x2 Pauli matrices. Using this Dirac algebra representation we see clearly the chiral decomposition 

^DD^ 
to2 + D^D 



where we have defined 



Thus, we write 



D = a^D^ , D^ = -a^D^ (5) 



r[A] = -^lndet(77i2 + Di:i^) - ^lndet(m2 + i:i^L>) 

= r(+)[A] + r(-)[A] (6) 

For later use, we recall the familiar properties of the matrices: 



,2 



m 



DD^ = -Dl + ]^F^,%,aa (7) 



D^D = -Dl + -F^,inl^aa (8) 

where 77°^ and 77°^ are the 't Hooft tensors O Unj . 

Now we make the following simple observation, that we can write the renormalized effective action as 

r,en [A] = 2r(±) [A] T (r w [A] - t[-J [A]) (9) 

Furthermore, we know that the difference of the renormalized effective action for the two chiralities takes a special 
form, related to the chiral anomaly: 

AFrenl^] ^ {v\tl[A]-TU[A] 

= ^(4^1^^)/^'^^-^- 

Thus, to evaluate the spinor effective action it is sufficient to evaluate the effective action for just one of the chiralities. 
That is, we can compute either F(+) or F^"^ but we do not need to compute both [T7] . 

This is a useful observation because there exist background fields for which the computation is significantly easier 
for one chirality than for the other. For example, if the background field is self-dual then F^^^fi'^^ = 0, since 77°^^ is 
anti-self-dual. Therefore the positive chirality operator reduces to the scalar Klein-Gordon operator: 

+ DD^ = w? -Dl (11) 

which implies that 

rspinor[^] = -2F,ealar[^] " AF[A] (12) 



This familiar result enables the computation of the quark determinant in an instanton background via the associated 
scalar determinant, which has a partial- wave decomposition [7l l8l [T2l [Ml [T5] . 
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In this paper we study a special class of background gauge fields that are not self-dual, but for which the Dirac 
operator still admits a partial wave decomposition [18H23]. Before coming to the radial decomposition, we first note 
that these gauge fields admit a simple chiral decomposition. Specifically, we consider gauge fields of the form 



(13) 



where rf^^, are the 't Hooft symbols FT^, and g{r) is a radial profile function, to be specified below. This type of 
background field is symmetric under 0(2) x 0(3) transformations, leading to a partial wave decomposition pOl I21j . 
This decomposition applies for any radial profile function g{r), so we have the freedom to study a wide class of 
background fields. In particular, we can investigate the role of zero modes, the presence or absence of which depends 
on the form of g{r). The associated field strength tensor is 



(14) 



Note that this field strength is not self-dual, due to the presence of the term proportional to g'{r). Thus, the positive 
chirality operator DD^ does not simplify as in (111. However, for this field the negative chirality operator does take 
a particularly simple form: 



(4g(r) +r.g'(r))a3 



(15) 



This is diagonal in spinor degrees of freedom, so if we work in the negative chirality sector, then we can immediately 
use our previous results for scalar Klein-Gordon operators [THU], just by including an extra "potential" equal to 
± (4(7(r) + r g'{r)). Thus, we can write 



D^D = 



4/ + 3 



dr - r'^g{rf - 4(7(r)Z3 - =f (4(7(r) + rg'ir)) 



= rn? 



n 



(16) 



where the quantum number I takes half-integer values: Z = 0, |, 1, |, 
In this way, we compute r'^~)[A], and then the full spinor effective action can then be obtained using (|9| and (10) 



while Z3 ranges from —I to /, in integer steps. 



III. THE PARTIAL- WAVE CUTOFF METHOD 



In this section we briefly recall the partial-wave cutoff method developed previously for radially separable fields of 
the form (13), adapted to the negative chirality sector of the spinor theory using (15). The basic idea of the partial- 



wave cutoff method involves separating the sum over the the quantum number I into a low partial-wave contribution, 
each term of which is computed using the (numerical) Gelfand-Yaglom method, and a high partial-wave contribution, 
whose sum is computed analytically using a WKB expansion. The regularization and renormalization procedure tells 
us how to combine these two contributions to yield the finite and renormalized effective action [Tllllj. 



A. Low partial-wave contribution 



The low partial-wave contribution for our system is given by 



r 



(-) 




det(m2 +H(^ij^^s)) 



det(m2 



-T / f rcc 



(17) 



where L is an arbitrary angular momentum cutoff. The factor — (21 



sum comes from adding the contributions of each spinor component in ( 16 ) 



1) is the degeneracy factor, and the s 
In order to evaluate this quantity, we 
use the Gel'fand-Yaglom method, which we now briefly review (for further details see [3{S]). Let Mi and A^2 denote 
two second-order radial differential operators on the interval r £ [0, 00) and let ^i{r) and 'E'2('") be solutions to the 
following initial value problem : 



Mi<^iir) = 0; ^.(r) - r 



21 



0. 



(18) 
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Then the ratio of the determinants is given by 



detA^, - hm r|4i ) . (19) 



det7W2 -R^oo \^2{R) 

In the negative chirahty sector we can take, Aii = m? + 'H(;./3 ,,.), and M.2 to be the corresponding free operator: i.e, 
the same operator with the background field set to zero: g{r) = 0. Thus, for a given radial profile function g(r) in 
( 13 1, for each value of [I, I3, s) we need to solve 

$± (r) + ^'^^'i (r) - (to2 + 4/3g(r) + r'^g{r f T [4.g(r) + rg'(r)]) $± (r) = , (20) 



with the initial value boundary condition in (18). The value of $ at r = 00 gives us the value of the determinant for 
that partial wave. In fact, the corresponding free equation is analytically soluble. It is numerically more convenient 
to define 

and solve numerically the corresponding initial value problem for S{r), as explained in [51 llOj. Then the contribution 
of the low-angular-momentum partial-waves to the effective action is 

r(-) = - f^W E [^+''^M + ^i''^^M]. (22) 

;=o,i,i,... '3=-' 

While each term in the sum over I is finite and simple to compute, the sum over I is divergent as L — >■ 00. In fact, 
only after adding the renormalized contribution of the high partial- wave modes do we obtain a finite and renormalized 
result for the effective action. 



B. High partial- wave contribution 

The high partial-wave contribution for our system is given by 



rlr^ = - E E ^(0 E 1- 



s=±i ;=L+i 



h=-i 



dct(m2 +'H(;^;3,s)) 

^det(m2+H[i-X,,))^ 



(23) 



Again, since by ( 15 ) we only need the negative chirality sector, and the negative chirality sector diagonalizes fully, we 
can apply our previous scalar analysis to each of the diagonal components, adding the appropriate term ±(^g + rg') 
to the Klein-Gordon operator. This modifies the detailed expressions as follows. In the large L limit, we have 



r 



(-) 



dr 



dr< 



8g{r)r^ 



f2 ^ 4 
3 



dr 



/2r3(3f3 + 8)g(r)^ 

(f2_,. 4)3/2 



6r^{5r^ + 28f^ + 32)g{rY 



45(f2+ 4)7/2 

-f 15(33^*^ + 335f'' + 1192f2 + mOO)g{r)^ 
+ 10r(15?~'^ + 184?'' + 776f2 + 1120) g (r) g' (r) 

+ 5r2(3f'^ + 38?-'' + 160f2 + 22A)g'{rf + 20r^{4: + P"^)"^ gir)g" (r)] 
r^{20g(rf + 10g{r)g'{r)r + g'{rfr'^) 
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r3(20.g(r)2 -|- lQg{r)g'{r)r + g'irfr^) 
\2e 



7 + 2 In L - 2 In 



2 + VWTl 



0\ 



L 



(24) 

(25) 

(26) 
(27) 

(28) 
(29) 
(30) 
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where f = ^ . We now identify the counterterm as 



<5r(-) = ^[i-7-21nM) / drr\20g{ry + 10g{r)g'{r)r+g'{rfr^), 



Where fj, is the renormahzation scale. Note that 

-J^ drr\F^,F^^) = dr r\%g{rf + Ag{r)g' {r)r + g' {rfr'') , 



1 poo rOO 

- / drr^iF^^F^,) = / dr r^8g{rf + Ag{r)g'{r)r) 
^ Jo Jo 

and thus, the counterterm corresponds to the following combination: 



^r^-) = ^ ( J - 7 - 21n/x j / dr r^U^.F^, + ^F^.^^ J . 



(31) 

(32) 
(33) 

(34) 



The appearance of the F^j^F^^, term is a special feature of the spinor calculation that does not occur in the scalar 
case. 

Having identified the counter-term, we can now write an explicit expression for the large L behavior of the high 
partial-wave contribution to the renormalized efi^ective action: 



V= / giog(r)lnL + ^g„(r)L" + ^g_„(r)— +0( — ), 

■''^ \ n=0 71=1 / 



(35) 



with the following expansion coefficients: 

Sg{r)r^ 



Q2{r) 
Qi{r) 
Q\os{r) 
Qo{r) 



3Vf2T4 
2r3(3f3 + 8)5(r)2 



1 



(f2 + 4)3/2 

r3(20g(r)2 + I0g{r)g'{r)r + g'irfr^) 



6r^(5f^-F28f^-F32)5(r) 



45(f 2 + 4)V2 



+ 15(33f^ + 335f'^ + 1192f2 + I600)g{r)'^ 
+ 10r(15f6 + 1847^-* + 776r2 + 1120)g{r)g'{r) 

+ 5r^{3f^ + 38f'' + 160f2 + 22A)g'{rf + 20r^(4 + f^fg{r)g"{r) 



giog(r) In 



2-h\/^2T4 



6r*(P+4f^)5(r)^ 



+ 2r2(f6 + I6f4 _^ + 128)g'{ry 

+ (-4f^ + 89f^ + 1104^* + +3456f2 + 5l20)g{rY 

+ 16r{2f^ + 21f'^ + 92f2 + lQO)g{r)g' (r) 

- 4r^(f*^ + 8f* + 16f'^)g{r)g"{r) 



(36) 



Note that involves L'^, L and InL terms that diverge as i ^ oo, but these divergences exactly cancel those 

of the Fl contribution (coming from the numerical Gel'fand-Yaglom computation for the low partial- wave modes), 
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yielding a finite renormalized result for the efi'ective action. This method is essentially exact (up to numerical precision) 
and accurately provides the value of the effective action for any value of the mass: 

r(.-„)(m)=4-)^™'^(m)+ri-)(m) . (37) 

The effective action calculated as above is finite for any non-zero value of the mass, however, from eq. (36), we note 
that Qo{r) contains a term proportional to In/i and therefore when we use "on-shell" renormalization {p — m), we 
have 



r'™(m) ^ y- J drQiog{r)j Inm ; m ^ . (38) 
In order to analyse the mass dependence of the effective action, we introduce a modified effective action defined as 

f-"(m) = r-"(m) + ( /"drgiog(r) I In/i. (39) 



which is independent of the renormaliztion scale and which is finite at m = 0. 



IV. RESULTS FOR SPECIFIC BACKGROUND PROFILES 



Within the class of background gauge fields (13), we still have the freedom to specify the radial profile function 
g{r). Here we note that the large r behavior of gjr) determines the presence or absence of zero- modes. To be specific, 
the integral of F^^F^i, counts the number of zero-modes. From (|33| we have: 

1 



2 



Jo 



drr^{F^,F^,) = / dr r^8g{rf + 4.g{r)g' {r)r) 



OO 

= 2{g{ryr ^ ■ (40) 

Therefore, as long as g{r) falls faster than we don't have zero-modes. In this paper, we study the background 
fields using two different profile functions to set the background field. The first one is 

gi(r) = B(l-Tanh(/3(VBr-0)), (41) 

where B, /? and ^ are parameters that control the amplitude, range and steepness of the potential. We may call gi{r) 
the "step potential". Since gi{r) falls off exponentially fast as r — >■ c», this case is free of zero-modes. 
The second profile function we use is 

ff2(r)^j.e-"'-V(p2 + r2)^ (42) 

where p and a are parameters that control the amplitude, range and steepness of the potential. This potential has 
the following properties: 



a > ^ / d'^xF^.^F^^ < oo , / d'^xF^.^F^^ = 



^4 p p 



a = =^ / d^xFn^Ff,^ oo , < 



< oo . (43) 



Note that g2{r) goes like l/r^ when we set a = 0, so that in this case there are zero modes. However, whenever the 
profile function g{r) decays as 1/r^ we have 



"'0 



drr^F^^F^, - / dr r^8g{r f + Ag{r)g' {r)r + g\r)\^) , 



^0 



oo 

drr'F^,F^,+ I drr\g'{rfr^) (44) 



which diverges logartihmically. In Figure [T] we show some plots of our chosen profile functions 51(7") and (?2('') for 
different values of their parameters. 

Our calculational method allows one to compute the effective action for arbitrary values of the fermion mass m. 
This provides us with a unique opportunity to study the validity of various approximate methods that yield estimates 
in the limits of large or small mass. In particular, we are able to probe exactly the to limit, which is of great 
interest in spinor theories but which is notoriously inaccessible by approximate means. 



FIG. 1: The upper figure shows the radial profile function: gi{r) = B{1 — Tanh(/3(v'Br — ^))). We have fixed B = 1 and 
P — 1 in all the curves. The parameter ^ which controls the range of the potential is varied. The lowest curve corresponds to 
5 = 1 (solid/red), followed by curves corresponding to ^ = 3/2 (dashes/blue), ^ = 2 (dots/black), ^ = 3 (dot-dashes/green). 
The lower figure shows a second radial profile function: 32 (r) = ue~"'~ /{p^ + r^). We have fixed u — 1 and p = 1 in all the 
curves, the lowest curve corresponds to a = 1 (solid/red), followed by curves corresponding to q = 1/20 (dashes/blue), and 
a = 1/400 (short-dashes/black). 



The large-mass expansion is the most general approximation method in the sense that that it may be apphed to 
calculate the effective action for any well-behaved background. Its main limitation is that it only applies for large 
values of the mass and, in fact, it diverges as m — >■ 0. Thus, it is not directly useful for probing issues related to 
massless quarks. To outline this method, consider for instance, the spinor case 



V. 



APPROXIMATION METHODS 



A. The large-mass expansion 
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In order to analyse the large-mass limit, we may take the small s limit in the proper-time integral and expand the 
trace factor in powers of s, yielding the heat kernel or Seeley-De Witt expansion: 



■■ Tre-(-^^)^^^^E^"«»[^M''] 



(45) 



The Seeley-DeWitt coefficients a„ are given by Lorentz traces of powers of F^^and its derivatives. For the Afj_{x) 
backgrounds given by (41) and (42) the Seeley-DeWitt coefficients are proportional to those found in the scalar case: 



0.1 
0-2 

03 



0, 
2 



^4 p p 



3 

"45 



(46) 



Since we are considering potentials of the form (131, with field strength ( 14 ) , we find an expansion of the renormalized 
effective action as: 



LM^2 



(47) 



The first two coefficients are: 
f(o) 

f 



LM 

(2) 



LM 



1 

6 7o 

1 

180 7o 



dr r^iSgirf + 4g{r)g'{r)r + g'irfr^) , 

00 

dr[24:r^g{r){l5g'{r) +r{9g"{r) +rg^^'>{r))) 



+r^{221g'irf + 9r^g"{rf + 2rg' {r){7lg" (r) + Grg^^^r)))] 



(48) 



Figure [2] presents a comparison between the effective action, as calculated using large-mass expansion and its exact 
value obtained using our partial-wave cutoff method. As expected, both methods agree well for large masses. 



0.5 



Eh 
X 
O 

IP-I 
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-2 
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FIG. 2: The graph shows the one-loop effective action, r(m), as a function of fermion mass m, with the radial profile function; 
gi{r) — B{1 — Tanh(/3(\/i3r — ^))). We have set B = 1, /3 = 1 and ^ = 1. The dots correspond to our exact numerical results, 
based on the Gelfand-Yaglom theorem, rGYT('Ti), and the solid line shows the outcome of the larg-mass expansion Ti^M{m). 
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B. The derivative expansion 



Another widely used approximation in gauge theories is the derivative expansion. This method is based on the 
fundamental result that for background gauge fields with constant field strength F^^^ it is possible to compute the 
renormalized effective action in a simple analytic form [^{5]. One can then expand around this soluble case, leading 
to an expansion of the form 



SoS - So[F] + S2[F, [dFf] + S4F, {dFf, [dFf] + ■■■ (49) 

The leading term in this expansion is the well-known one-loop effective action for constant backgrounds, first computed 
by Heisenberg and Euler [3 . In euclidean QED, the corresponding one- loop effective Lagrangians for spinor and scalar 
QED are 

£spi„or(a,&) - Je-"'"{a5s2coth(as)coth(6s) + l- y (a2 + &2)|, (50) 

and 

where ztia and ±6 are the eigenvalues of the 4x4 constant matrix F^^, and are related to the field invariants in the 
following way: 



a' + 6' = ^F^.i^"'' , ah=]^F^,F^^'' . (52) 



The leading term in the derivative expansion is obtained by simply replacing the the constants a and h by the 
corresponding space dependent quantities inside these integral expressions, and then obtaining the effective action 
from the effective Lagrangian by integrating over spacetime. For the class of backgrounds that we are considering, 
the corresponding substitution is 

a — > a{r) = 2g{r) , 

b — > b{r) = 2g{r) + rg'{r) , (53) 

such that, 

POO 

rDE-2^M r''drCiair),bir)). (54) 

JQ 

At first sight, one would expect the derivative expansion to be a good approximation at large mass, similar to the 



large mass expansion, since the expansion over derivatives in (491 is balanced by inverse powers of m. However, 
the situation is more subtle than that naive expectation, as the derivative expansion expression at a given order is 
a resummation of all powers of the field strength with a fixed number of derivatives. Thus, one might expect that 
the derivative expansion is better than the large mass expansion for smaller values of the mass. Indeed, in previous 
work on the partial-wave cutoff method in scalar theories it was found that even the leading term of the derivative 
expansion provides a surprisingly accurate approximation to the mass dependence of effective action |10j . even for 
very small mass. In the next section we study this question for spinor theories and find an interesting difference. 



VI. APPROXIMATION METHODS VERSUS EXACT CALCULATION 



In this section we use exact results obtained from our partial- wave cutoff method to probe the validity of the large- 
mass and deriative expansion approximations, for spinor QED. In order to emphasize the similarities and differences 
between the scalar and spinor cases, we also show some results for the scalar action as well. 



A. The scalar case 



To exemplify how the partial-wave cutoff method works in the scalar theory, we use the background field given by 
(7i(r), choosing different values of the range parameter ^. We present a comparison between the the exact effective 
action as calculated with the partial-wave cutoff method, the large-mass expansion and the derivative expansion. 
From the plots in Figure |3] we make the following observations: 
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• The large-mass expansion result agrees very well with the exact result, already for m 2. 

• The leading order derivative expansion is good at large to, and also provides a very good approximation to the 
effective action for large values of the parameter ^. 

• In particular, as long as the steepness parameter ^ is large enough, the derivative expansion provides accurate 
results in both the large-mass and small-mass regimes. 



B. The spinor case 



In this section we compare and analyse the different approximation methods against the GYT method. The first 
background configuration we investigate is the plots shown in fig. (|4| correspond to the GYT method, the 

large-mass expansion and the derivative expansion. From the plot in Figure|4j we make the following observations: 

• The partial-wave cutoff method produces a finite value for the effective action in the small-mass regime. 

• As in scalar QED, for spinor QED the large-mass expansion result agrees very well with the exact result, already 
for m ^ 2. 

• As in scalar QED, for spinor QED the leading order derivative expansion is good at large m. However, the 
leading order derivative expansion result diverges in the small-mass regime. 



To understand why the leading-order derivative expansion fails in the small-mass limit, we analyse equation (50 1 in 
the small-mass regime. The small- mass regime is obtained by taking s — > oo, which gives: 



spinor 



(a, 6) 



1 

1 

4^ 



ah-\{a' + b') 



ds 



ab - 





Inm 



(55) 



We recognize one In to divergence proportional to F^^F^^, and another proportional to F^^F^^. Contrast this with 
scalar QED where 



scalar 



(a, 6) 



1 



167r2 
1 rl 



Q+Ua' + b') 
b 



ds 



In TO 



(56) 



which has a Iuto divergence proportional to Ff^^F^^j, but none proportional to F^^Ffj_,y. This is simply a reflection 
of the fact that zero modes can occur in spinor QED but not in scalar QED, with the F^^Ffj_i, term determining the 
number of zero modes. 
But note that £, 



-'spinor V 



(a, 6) = >Cspinor(— a, ~^)- Therefore, (55 1 really should read 



-■spinor 



(a, 5) 



1 

4^ 



\ab\-l{a' + b') 



Ini 



(57) 



Therefore, our naive application of the leading-order derivative expansion includes a term in the effective action that 
behaves in the small mass limit as 



/ 



d'^x\a{r)b{r)\ 1 Iuto 



(58) 



instead of the correct form 



d'^xa{r)b{r)j\Tim (59) 

It is this latter form that counts the number of zero modes, and appears in the counter-term. To show the effect 
of this, consider the radial profile function gi{r) for which J d^x F^^^F^^ = 0, indicating the absence of zero modes. 
For this profile function, this integral vanishes because the integrand changes sign. But the derivative expansion 

expression does not allow for such changes of sign, and so we effectively compute / d'^x Ffj^^F^^ , which is non-zero. 
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This mis-match is demonstrated clearly in Figures [5] which show that the divergent small-mass behavior of the 
derivative expansion corresponds exactly to 

foE - /M = J Hr)h{r)\^ Inm ; m ^ . (60) 

The number of zero-modes is given by J d'^xa{r)b{r) and it is always equal to zero for the backgrounds of the class 
gi{r), however, ^ J d'^ x\a{r)b{r)\ does not vanish and this adds an incorrect residual logarithmic dependence to the 
derivative expansion, as we have shown. No such mis-match occurs for scalar QED because the small m behavior of 
the derivative expansion expression only involves the combination i^^^-F^i, oc [a^(r) -|- b'^{r)], which is always positive. 
Thus there is no Inm divergence in the small m behavior of the derivative expansion plots shown in Figure [3] for scalar 
QED. 

The second background wc examined is the one given by the profile function g2{'f)- Setting different values of the 
parameter a, we can control the rate of decay of the potential. In figure [6] we present a comparison between the 
different calculation methods for this kind of background. We corroborate once more how the derivative expansion 
fails at the zero-mass limit. Also, the derivative expansion shows better accuracy in a wider mass-range as the for 
those fields with slower variation (small a) , as expected. It also evident that the residual logarithm is not the dominant 
term in the derivative expansion when we move from the small-mass regime in to the large-mass regime. 

VII. CONCLUSION 

We have extended the partial-wave cutoff method to spinor theories with nontrivial [and non-self-dual] radially 
symmetric gauge backgrounds. Different background fields were tested, resulting always in accurate values of the 
effective action for both the large-mass and small-mass regimes. We provided an example of how our method allows 
one to systematically investigate how the effective action responds to different characteristics of the background fields, 
such as range, amplitude or rate of variation. We have also analyzed how certain approximation methods compare with 
these exact results, in the different mass regimes, also comparing the scalar case with the the spinor case. Specifically, 
we have tested the large-mass expansion and the derivative expansion. We have shown that the large mass expansion 
works extremely well, as in the scalar case. However, the derivative expansion behaves in a different manner in 
the small-mass regime for the spinor theory. We have explained this fact, both qualitatively and quantitatively, as 
resulting from the changing sign of the local quantity Ffj^nF^u, and which in the derivative expansion approximation is 
assumed to be constant and therefore of fixed sign. We have also shown that this is directly related to the appearance 
of fcrmion zero modes. 

GD and AH were supported by the US DOE grant DE-FG02-92ER40716, and HM was supported by the National 
Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (No 2010- 
0011223). 
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FIG. 3: The graphs show, for the case of scalar QED, our exact effective action, fGYT(?Ti) (big-dots/red), the derivative 
expansion expression TnEim) (small-dots/blue), and the large mass expansion expression, rLM(m-) (solid-line). We use the 
profile function gi{r) — B{1 — Tanh(/3(\/B'' — C))) with B = 1, /3 = 1. We show results for three values of the range parameter: 
^ = 1 (upper graph), ^ = 3/2 (middle graph) and 5 = 2 (lower graph). Notice that the derivative expansion is a reasonable 
approximation even at zero mass. 
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FIG. 4: For spinor QED, this plot shows the exact effective action, f GYT(?n) (big-dots/red), the derivative expansion expression, 
TuE{m) (small-dots/blue), and the large mass expansion expression, ^LM(n^) (solid-line). We use the same profile function as 
in the upper plot of Figure [s] [which is for scalar QED]: gi{r) = B{1 - Tanh(/3(^r - (,))) with B = 1, /3 = 1 and 5 = 1. Note 
the very different behavior at small mass compared to the scalar QED case. 
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FIG. 5: For spinor QED, these graphs show the exact effective action, fGYT(fTi) (big-dots/red), the derivative expansion 
expression, rDE("i) (small-dots/blue), and the large-mass expansion expression rLM(»Ti) (sohd-hne). The dashed hne represents 
the residual logarithm f{m) = d*x\a{r)b{r)\)\nm. We use the profile function gi(r) — B{1 — Tanh(/3(\/i3f — C))) with 

B = 1, P — 1. We show results for three values of the range parameter: ^ = 1 (upper graph), ^ = 3/2 (middle graph) and C = 2 
(lower graph). Note that the divergent behavior at small mass is well fitted by the residual logarithm, as described in the text. 
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FIG. 6: For spinor QED, these graphs show the exact effective action, TGYTim) (big-dots/red), the derivative expansion ex- 
pression, TuEim) (smaU-dots/blue), and the large-mass expansion expression rLM(w.) (sohd-hne). The dashed hne corresponds 
to ^{J d'^x\a{r)b{r)\)\nm. We use the profile function g2{r) — ve j (p^ -|- r^) with v = \ and p —\. We show three values 
of the decay rate parameter: a = 1 (upper graph), a — 1/20 (middle graph) and a — 1/400 (lower graph). Note that the 
divergent behavior at small mass is well fitted by the residual logarithm, as described in the text. 



